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OPTIMAL DESIGN OP WATER DISTRIBUTION SISTM 

R.S. SOLANKI 

Em'iroriiiiental Engineering Division 
Department of Civil Engineering 
Indian Institute of Technology, Eanpur 
-July., 1980 

MpstLof-rthacmatiaods developed so far for the optimal 
design of looped water distribution systems employ nonlinear 
optimization techniques along with a network solver. 
E.onlinear optimization techniques require a high degree of 
experience in applying these to complex problems in order to 
give initial values of various parameters. It becomes very 
difficult to solve the nonlinear problem when measures are 
taken to account for the variable demand pattern and discrete 
nature of some of the variables (e.g. pipe diameters). 

The present work is aimed at developing a computer 
program based on the linear Programming Gradient (IPG) method 
as proposed by Alperovits and Shamir, . No network solver is 
required in this approach. The progi^am develop©3 here deals 
with the problem of designing an optimal water distribution 
system for a single loading pattern and gives the solution 
in terms of the available pipe diameters, head at the 
source (elevation of the reservoir or the horse-power of the 
pump), and capacity of the booster pump if required at a 
specified location. 


1 . IHTBOPUCTIOH- 


Water is a primary necessity of life. Throughout 
recorded history large cities have been concerned with their 
water supplies. Even ancient cities found that local 
sources of supply - shallow wells, springs, and brooks - 
were inadequate to meet the very modest sanitary demands, 
and the inhabitants were constrained to build aqueducts 
which could bring water from distant sources. In the present 
times, the need of supplying treated water to community is 
well recognized. Extensive distribution systems are 
needed to deliver water to the individual consumers, 

A water distribution system is meant for 
supplying water at sufficient pressure and inadequate 
quantity to various points in a locality. The distribution 
system consists of various components such as pipes, 
reservoirs, pumps, and valves. An engineer is concerned 
with the design, operation and maintenance of a distribution 
system i.e. deciding about the sizes or capacities of 
different components in the systm and their operating 
policies. 

There are some basic rules forming ^ 

foundation of all design procedures for a water distribution 
network, such as flow continuity at the nodes, head conti- 
nuity along paths of the network, and the head loss vs. 
flow relationships for different components in the system. 
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In. the conventional design procedures, tentative sizes 
of components in the system are ^sumed and then the 
analysis techniques are applied to check the feasibility 
of the system regarding the required amount of water being 
supplied at adequate pressures to all nodes of the network. 
Hence, a number of combinations of the component sizes 
could be possible to serve the same purpose. Optimization 
techniques are very useful in selecting the best alter- 
native out of a number of possible solutions in view of the 
constraints on the resources and the objective to be 
achieved. 

Considerable saving in cost can be achieved by 
using the optimization techniques over the conventional 
methods of water distribution system design and operation. 
This water distribution system is often the major invest- 
ment of a municipal waterworks. Such saving in cost is, 
therefore, of immediate concern as India being one of the 
participating countries at the World Water Conference 
held in 1977 at Mar del Plata, Argentina, has committed 
itself to the drinking water supply and sanitairion programme 
during the decade 1981-1990, Ho citizen of this country 
should be denied access to safe potable water by the year 
1990. This would involve investments well over Rs. 15*000 
erores . Development of efficient techniques of design 
could go a long way towards the fulfilment of a programme 
of such a vast magnitude . 



1.1 Objective of the Study 


Various investigators have developed a number of 
optimization algorithms for designing a water distribution 
system. Most of methods proposed utilize nonlinear optimi- 
zation techniques for looped network systems. However, 
recently some work has been done in the area of application 

of linear optimization techniques for designing a looped 
2 

network , Application of linear programming may be 
promising in view of its being more efficient compared to 
nonlinear programming. 

The present work is aimed at developing a 
computer programme based on 'linear Programming Gradient' 
method developed by Shamir and testing it on a few networks 
of varying sizes. The various cost data taken are in 
accordance with the local conditions. 
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2.- LITBRATURE REyiSW 

2 . 1 Mathematical Models for Steady State Analysis 

Water distribution systems are designed to 
deliver water from sources to consumers in the required 
quantity and under a satisfactory pressure through pipeline 
networks equipped with a variety of components such as 
pumps, valves of various types, distribution reservoirs. 

The mathematical model of a network consists of links 
connected at nodes. Nodes are points at which two elements 
connect, or at which a flow enters or leaves the network. 

To maintain mathematical tractability it is customary to 
model only the major structure using some appropriate 
representation of the detailed structure of the actual 
system. Bor example, it is customary to include in the 
model only pipes above some minimum diameters, and also to 
lump withdrawls taken along the pipes and assign them to 
the nodes. 

At each node we define the head and the consumption. 
Head is the sum of the topographic elevation, the pressure 
and velocity heads (the latter is usually negligible) 

While there exists a head at any point along a pipe we 
shall refer by "heads” to the set of values at the nodes. 

At certain nodes there are imputs of water into the system. 
At others there are cor^umptions , which are viewed as 
negative inputs. 
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Bach link is chacacierized by a physical law 
which relates the flow through it to the head difference 
hetween its ends. Bor pipes there are several empirical 
flow equations; one commonly used is the Hazen— Williams 
equation I 

where Q is the discharge; smoothness coefficient 

(the Hazen-Williams coefficient); 3) the diameter; the 
head difference between the two end nodes, and 1 the length, 
o<; is a numerical coefficient, whose value depends on the 
units used (For Q in m^/sec and D in cm < = 5.4 x 10 
for Q in cfs and 3) in inches ci = 6,28 x lO"'^), £^H/1 is 
the hydraulic gradient (hereafter denoted by J). 

Bor centrifugal pumps, the head added by the pump, 
is usually approximated by a pol;^omial of the form 

= a + bQ + cQ^ (2) 

Every other type of element has its own law. Reservoirs 
are connected at nodes; water level in the reservoir is 
the head at that node. 

The mathematical model of the network is a set of 
simultaneous non-linear algebraic equations, which corres- 
pond to a steady state flow in the network under fixed 
boundary conditions? a set of inflows and consumptions at 
the nodes and a set of fixed beads at specified nodes. 

The solution of these equations is called a steady state 
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flow solution. There are two basic types of equations? 
node equations and path equations. The first express 
material continuity at a node. For node 3? 

where is the flow from node i to node 3 and I . is the 

inflow into node 3. Consumptions (withdrawls) axe denoted 

by C . and appear as negative values of I . . 

'' 3 

Path equations equate the head difference, b^, 
between the end nodes of path (p) in the network to the sum 
of head gains and losses in all links belonging to this path? 




= b 


( 4 ) 


A path may connect any two nodes. Usually path equations 
are formulated between pairs of nodes at which the heads 
are loiown. A special case are loop equations, where the 
two ends are at the same node, and then b^ = 0. 

There are many different ways to construct the 
mathematical model. In a network with N nodes, a set of 
U node equations fully determines the flow solution, 
provided one head (a reference head) is given. When loop 
equations are used one needs as many equations as there 
are basic loops, i.e,, loops which do not have pipes inter- 
secting them. General path equations can be formulated in 
many v?ays and the rule is that they should be mutually 
independent. Any one of these formulations can be used to 
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obtain the flow solution i.e. the heads and flows throughout 
the network. The set of simultaneous non-linear algebraic 
equations can be solved by an appropriate iterative 

technique - Hardy-Cross, Newton-Raphson, or lineariza- 
tion^^ 

When there is internal storage in the system, i.e. 
there are operational reservoirs then the sequence of flow 
solutions for varying withdrawls over time (e.g. a day) is 
of importance. Subsequent flow solutions are linked 
through the changes in water levels in the reservoirs. 
Simultation of such a system's operation is carried out by 
obtaining a flow solution for the initial boundary condi- 
tions, the resulting flows are used to update reservoir 
levels, then these levels are used as boundary conditions 
for the next flow solution, and so on^*^^ 

2,2 Problems of Planning, Design and Operation 

The layout of the small-diameter pipes in a 
distribution network is essentially fixed by the land 
development. Water has to be delivered to every building 
and pipes usually follow streets, so the layout of the 
finer grid is fixed to a large extent. On the other hand , 
the engineer has to select the pipe materials and diameters, 
consider the possibility of breaking network into pressure 
zones (separated by special pressure regulating devices), 
allocate the loads to the various sources, deteimine the 
layout and sizing of the feeder mains, and fix the 


8 


locations and designs of the pmps, reservoirs and other 
facilities. Planning is the phase of selecting the layout 
and Eiain features of the system. Design is the phase of 
fixing sizes and characteristics of the various components. 
The two are complementary and should he carried out simul- 
taneously. 

Distribution systems have to operate under time- 
varying conditions. Considering long range changes, on 
the scale of years, there is the need to meet increasing 
demands . This capacity expansion problem is common to many 
areas of engineering. The scale of time-varying conditions 
which is of specific concern is associated with the 
operation of the system. Water distribution systems operate . 
under loads which change over the day, the week and the 
seasons. Sources, pumps, valves and reservoirs are operated 
to meet these varying demands , Setting the operating policy 
of the system is an integral part of the design process, 
since sizing the components depends on their operation 
and vice versa. 

In the design one would like to consider expli- 
citly the detailed operation of the system, that is, the 
hour by hour position of the pumps, valves and reservoirs. 

No way has been found as yet to formulate and solve the 
design and optimization problem in this form. Instead one 
can Inelude in the foimulation of the design problem speci- 
fication of the operation under one or several loadings 
(sets of consumptions) . Hiese are ’‘typical” or "critical" 
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loadings - for example, average daily consumptions, maximum, 
demands during the day, low demands (which usually occur at 
night, during which reservoirs can he filled), high demands 
for fire fighting (usually concentrated in the highly 
congested business district), etc. The solution of this 
design and operation problem provides both sizing of the 
various components and their operation under these typical 
loadings . 

Planning, design and operation are thus three 
aspects of a single problem, but under certain circumstances 
the task may be narrowed to design and operation, or 
operation alone. In the following section which deals with 
optimization, we shall treat planning and design as a single 
problem. The method is as follows: one specifies facilities 
(pipes, pumps, valves, reservoirs) wherever they seem 
reasonablej the optimization procedure is allowed to set 
any design variable to zero, thereby eliminating the 
element from the solution. This procedure will not create 
a facility where one was not specified, and is therefore 
limited to those configurations stipulated by the designer. 
Still, this procedure allows for selection among alternatives 

2.5 Formulation of Network Optimization Models 

We shall first discuss the general structure of 
optimization problems for planning design and operation 
of water distribution systems and then go on to describe 
and discuss the methods which have been developed to solve 


them. 
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2.3.1 Decision yarlaible a 

In. selecting pipes one has to decide on their 
material, diameter and wall thickness (pressure hearing 
capacity). In the optimization, it is usually asstimed that 
material (i.e. smoothness) and wall thickness have been 
fixed, and the remaining decision is on the diameter. In 
selecting pumps there are various considerations, but for 
the optimization only the head vs. discharge character- 
istics of each pump will be considered as the decision 
variable. Similar considerations bold for valves and 
reservoirs . 

Decision variables of the planning and design 
problem includes 

a) Pipe diameters, 

b) Pump locations and characteristics, 

c) Valve locations, and 

d) Reservoir locations and sizes. 

Vtoen operation of the system under a set of 
tjrpical or critical loadings is included in the formulation, 
the following decision variables are added s 

e) Pumps to be operated (ON/OPP) under each loading, and 

f ) Valves to be operated (amount of pressure loss provided 
by the valve) under each loading. 

2.5.2 Objective Punotions 

Minimum cost is the criterion most often used 
in optimization of water distribution systems. Total cost 
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is made of capital plus operating costs. Ihe latter 
usually reflects only energy costs, since operation and 
maintenance may be included in tbe capital cost. 

Performance indicators, such as minimum 
pressures at supply nodes, are usually treated as const- 
raints in the optimization. Some work has been done in 
which certain performance criteria were used in the objective 
functions in a multiple-objective formulation^^. 

2.5.5 Constraints 

Several types of constraints appear in the 
optimization models. First, the physical laws of flow in 
the network have to be satisfied. These are equality 
constraints for continuity of mass and/or of hydraulic 
head lines, equations ( 5 ) and/or (4). The consumptions are 
usually treated as fixed externally, so they appear as 
constraints in the node continuity equations. Limits are 
normally set for the heads or pressures at some nodes. 

Minimum pressures are to be guaranteed under all loadings, 
to meet service standards for domestic and industrial 
consumers and to ensure sufficient operating conditions 
for fire fighting. Maximum pressures are specified when 
there is a danger of pipes bursting or equipment being 
damaged under excessive pressure. 

Special types of constraints may arise from 
specific formulations of the optimization models. These 
will be discussed specif ioally for each model to be presented 
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2,4 Branching letwork Models 
2.4.’5 Linear Prograonming Model 


Considering a branching network, supplied from 
one or more sources by gravity, and designed for a single 
loading, having specified consumptions at the supply nodes 
to be satisfied ( 0 - at node j)? at some or all of the 

ti 

nodes the head, H., is to be within a given range, 

to HMAX.. The layout is given, and the length of the link 
3 

connecting nodes i and 3 is L^.. 

. 16,17,19,24 

The linear programming design procedure 

is based on a special selection of the decision variables . 

for each link allow a set of "candidate diameters , the 

decision variables being the lengths of the segments of 


these diameters within the link* Denoting by the 

length of the pipe segment of the m—th diameter in the 


link between nodes i and 3 , then 

T X. .„ = L.. for all (i, 3 ) 

m 

where each link may have a different set of candidate 
diameters. For a branching network in which the consump- 
tions are known, the discharges in all links, 
fixed . The head loss in segment m of the link is 




S j . . X. 

13 m 131 a 


for all Cifa.’a) 


( 6 ) 


Where J is the hydraulic gradient ( ^H/L in eq. (1)) which 
is a function of the discharge and of the diameter (if the 
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smoothness is assumed to be selected in advance, and there-- 
fore fixed). 

Starting from any node in the system, s, at which 
the head is fixed (e.g, a reservoir), and selecting any 
path from it to node n, at which the head has to be within 
a given range, one may formulate the constraint 


HMIN < 
n 




(7) 


The first summation is over all links along the 
selected path, and the second over all segments of the link. 
The signs of the terms depend on the direction of flow in 
the link. In order to reduce the number of constraints and 
improve computational efficiency, head constraints may be 
formulated only for part of the nodes, provided they suffice 
to ensure that pressures throughout the network are within 
their acceptable ranges. If this method is used one has to 
exaainc the solution to ascertain that all heads are 
satisfactory? wherever they are not, a new constraint has to 
be added and the problem resolved. 

The cost of a pipeline with a fixed diameter can 
reasonably be taken as linearly proportional to its length. 
Thus, the total cost of the pipeline network is 


(i»3) ® 





Minimization of (8), subject to (5), (7) and non- 
negativity of the X is a iinear program (LP). 
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It can be shown that if the cost of a pipe is a 
convex function of the diameter (as it noimially is) then 
in the optimal solution of the IP each link will contain at 
most two segments, their diameters being adjacent on the 
candidate list for that link^. 

The LP formulation can be extended to include the 
cost of pumps, and their operation, using linear or linearized 
cost functions * . Reservoirs can also be included, using 

the head in the reservoir as the decision variable, and 
fitting it with a linear cost function. More than one 
loading can be considered, Bach loading results in a set of 
constraints of type (7), possibly with different bounds on 
the heads for each loading, and the entire set is solved 
simultaneously in the LP, If energy costs are included the 
objective function contains a weighted sum of the energy 
costs for operating under the different loadings. 

2.4.2 Dynamic Programming Model 

Optimal design of a branching network, with or 
without pumps and reservoirs, can easily be formulated as 
a dynamic programming (DP) problem. The solution requires 
more computer time than the LP method, but the formulation 
is free of certain shortcomings present in the LP. 

pi . 

Rally used DP to optimize the diameter and 
wall thickness (which determines the pressure bearing 
capacity) of the segments of a pipeline, as well as the 
heads added by the pump located along it. The objective 
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function included capital pltjs operating costs. Similarly, 
Liang used DP to optimize the diameters of segments 
between takeoff points to consumers along a pipeline fed "by 
a pump at Its upstream end. The objective was to minimize 
capital cost, subject to minimum pressure constraints at 
the takeoffs . 

Probably because of its computational inferiority, 
the use of DP for branching networks has not been developed 
to an operational stage. Still, for completeness, the 
structure of DP formulation is presented here. 

The network is divided into segments, A segment 
is defined between adjacent takeoff points, so that the 
discharge remains the same along the segment, or, if these 
segments are too long, a finer division may be used, 
allowing pipe properties to change from one segment to the 
next.- Takeoffs are given, and minimum pressures are to be 
satisfied at each takeoff node. Decision variables may 
include the diameters, material (roughness and strength) 
and class (pressure bearing capacity) of pipes, and the 
capacities of pumps, for clarity of the presentation we 
shall assume that material and class of the pipes have been 
fixed, so that diameters are the only decision variables. 

The objective function may include capital cost of pipes 
and pumps, energy costs, and any benefits, costs or penalties 
which are the functions of the heads at the nodes. No 
assumptions, such as continuity or comrexity, have to ^ 
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made about these functions. The quantities to be supplied 
at nodes are assumed fixed, and constraints may be imposed 
on the minimum and/or maximum beads at the nodes. 

The state variables are the heads at the nodes, 
and the nodes are the stages. Computation proceeds upstream, 
starting from the downstream end of each branch. The 
recursive equation of the DP is : 

= ““ L®‘V * 

Djj. is the diameter of segment k which connects node ( j+1 ) 

to its downstream neighbour, node 3. g(Dj^) is the cost of 

this segment. is the head at node (J+ 1 )j H. is the 

O/vvCt 

head at node 3 given and the diameter can easily 

be computed since the discharge in the segment is known 
f(H.^-, H . ) is the cost (or benefit) associated with the 
link, for the given heads at its two ends. P^(H.) is the 
optimal value for the poarfcion of the system downstream of 
node 3, given the head at node 3. The minimization is over 
all admissible values of Dj^., and is performed for each of a 
set of discretized values of over its admissible range. 

Whenever , Djj., C^) is outside its admissible 

range, the examined D^j. is disallowed. 

Por a pimip located between nodes and 3 

one /uses : ■ \ \ ' 


,) = Min [gCH.*,, Hj) * Hj) * ( 10 ) 


H. L 
3 
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+ here is the capital plus operating cost for a 

pump designed to deliver the l£nov?n discharge, Qj^., fr<aa 
head at its intake to head H.( ) at its discharge. 

The minimization is over a set of discrete values of H. . 

tJ 

At every branching node of the network one adds 
# 

up the values of for the downstream hramhes which 
connect at it. At such a node j, pt(H ) is still the 
optimal cost of the part of the network dovixstream from it, 
except that now it is a sum of the optimal costs for all 
branches originating at node Thus, one follows the 
single line procedure outlined above, starting from all 
downstream extremities of the network. 


2.5 Looped Network Models 

2.5.1 Linear Programming Models 

26 

Lai and Schaake have attempted this problem by 
making the assumption that the beads at all nodes, as well 
as the demaMs are given in advance. The solution thus 
gives the optimal diameters for the assumed pressure pattern. 
Since the beads are fixed, the flow through each link is a 
function of the link properties only, and if the length and 
pipe material are given, the flow is a function only of the 

diameter. Por node j at which a demand G. has to be satisfied, 

. 4 : 

the following constraint is a statement of continuity at 
the node ; 


*1 

i 


% 

i 


ii 13 


s G. 
3 


(11) 
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where is a ooefficien.t whose value is determined bj the 
given datas the heads at nodes i and ;j» and the length and 
smoothness of the pipe connecting them. The objective 


function Lai and Schaake used was 


51 

(i, j) 


a I. . 
ID 


ID 


* b 


f 


(i»D) 





( 12 ) 


a and b are constants which account for unit conversion » a 
present value factor, etc. e is a coefficient whose value 
is determined through analysis of pipeline cost data. The 
first sum is the capital cost of the pipelines, the second 
is the cost of energy lost in flow through the pipes , and 
the third is the cost of the residual energy at the supply 
nodes. There are no pumps within the network, and the 
energy terms represent the cost of supplying water to the 
network from the external source by a pump. By makir^ the 
substitution Y.-.- = the constraints (11) become linear. 

Using eq. (1) for the flows, the objective function (12) 
becomes s 



(i»D) 


(a L. . + b. - Y. .) 

^ 13 13 13 ^13' 


(13) 


where the coefficients b^ ^ are based on the heads and flows. 

Minimization of (13) subject to constraints of type (11) 

for all supply nodes and to non-negativity of Y^^ is 

linear program (which Lai and Schaake solved by a 

developed iterative LP program). The method was used in a 

10 . 

study of New York City’s primary distribution system in 
which several performance criteria v^re also introduced as 



objectives (e.g, sxm of tbe residual pressures at all the 
supply nodes, the residual pressure at the farthest supply 
point, etc.). 

Kally approached the problem by extending the 
IP formulation for branching systems. His reasoning is as 
follows. In a branching system, if one changes the length 
of a pipe segment (which has some fixed diameter) the 
resulting changes in the heaSs at nodes are linearly 
proportional to the magnitude of this change. In a looped 
network the same effect is non-linear, due to the redistri- 
bution of flows once the design is changed. Still, if the 
change in length of a segment is small enough the resulting 
change in heads is approximately linearly proportional to 
the change in length. The ultimate decision variables in 
Kally 's fomulation are the lengths of the segments, the 
same as was for a branching network, and the objective 
function is also the same. The problem is solved through 
a sequence of linear programmes 5 the decision variables in 
each LP are the lengths in each link along which the diameter 
is to be changed from one diameter to another, i.e . the 
length to be taken away from one segment and given to 
another. Denoting by the length in link (i,j) of 

change from the present diameter to diameter m (of a candi- 
date list), these changes have to satisfy 
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that is, the sum of all changes cannot exceed the total 
length of the link. Also, considering the minimum heads 
required at certain nodes, HMIN. , the changes in diameters 
are limited by 


(i,3) m 






( 15 ) 


where is the head at node k in the present iteration, and 
(Hj^ - is therefore the extra head which can still he 

eliminated (if feasible from other considerations). The 
coefficients (^ are the linear approximations for 

the rate of change of head with respect to changes in segment 
lengths. Kally’s approach was to obtain these values as the 
difference in heads at the nodes between two solutions - one 
with existing segments, and another with one = 1, 

For each As - = 1 one has to run a network solver and 

obtain a flow solution, then one computes the head differ- 
ences jii^hj^ between it and the "basic” flow solution - the 
one with ^ s^^^^ =0* 

The objective function for each iteration's IP is 

Min ^ ^ C-.s.. (16) 

i7o ~m 

where C- - is the cost of changing one unit length of pipe 
in link (i,j) from its present diameter to diameter m. 

Upon solution of the IP, the changes of diameters 
over lengths introduced and a new iteration is 

begun: a "basic" flow solution is computed, a series of flow 
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solutions for all A = 1 are obtained and (16) is mini- 
mized, subject to (14), (15) and non-negativity of the s, 

Kohlhass and Mattern^^ used linear programming in 
optimizing a looped network in which the heads are fixed in 
advance - a condition similar to the one stipulated by Lai 
and Schaake^^. 

Recently, a more general method has been developed^ 
and introduced into practice. The method called Linear 
Programming Gradient method (LPG) is based on the following 
reasoning. If the flows throughout a looped network are 
known then its optimal design can be obtained by an LP 
fomulation similar to that for a branching network. Ibe 
optimization has, therefore, to find the optimal flow 
distribution, and for it the optimal design. This is achieved 
by a hierarchical approach. In the lower level of the 
hierarchy the optimal design for a particular flow distri- 
bution is obtained by LP; in the higher level, the flow 
distribution is modified, using certain results of the LP 
solution, towards an optimal flow distribution. This 
procedure is continued iteratively until some termination 
criterion is met. 

The LP for a fixed flow distribution in a looped 
network operating under gravity for a single loading iss 
minimize (8),subject to (5) and (7), to non-negativity of 
the X, and to the additional constraints 


2 - 

(i,j)€p 


z. 


^ijm ^ijm 



(17) 
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where p designates a path in the network, and 13^ is the 
(known) head difference between its ends (all other notation 
is the same as in Section 4.1). Eq. (17) has to hold for 
all loops in the network with b = 0. When pumps, valves 

Jcr 

or reservoirs are to be included, the objective function and 

constraints have to be augmented, as will be explained below. 

Having added constraints (17), the IP can be solved, 

and the set of optimal segments will be such that the 

network is hydraulically balanced by virtue of the fact 

that the constraints (17) have been satisfied. Denoting by 

Q the vector of flows in all the paths, which may be any 

arbitrary set of flows as long as they satisfy continuity 

at all nodes, then the optimal cost of the network, P, for 
— % ■ 

this Q can be written as 


P = LP (t. ) 


(18) 


where IP denotes that P is the outcome of a linear program. 
Next Q is modified in a way which approaches optimality. 
Denoting by the change of flow in path p, then 

ir 

2 p 




DP ^ “p ^ 3 P 3 ^r 

r€R 


= W. 






P * 'r fcH 


b 


(19) 


where Wp and \4j, are the dual variables of constraint (17) for 
the paths p and r, respectively, where R are the paths which 
share a link (or more than one link) with path p. Usi)^ 

(17) and the definition of J from eq. (1 ): 
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^ 2. 1.852 o9-.852 571 .852 

{i.d)^;P n ^3“ 


DTf -S'? X. , 

^3^ 13221 


’•®52 ,.^- jr- ^ AH 


ionx 


( 20 ) 


9(ZSQp)= ^(0p) because both are increExental changes in 
the flow in the path. ard bave already been used 

in setting up the LP, so once it has been solved and the 
duals Wp and are known, the components of the gradient. 

r 1 


J P 

a^^Qp) 


= 1 .852 






idm 


S' 


w„ 


r£.R (i,d)6.r m 


Q-- ^-^ijm 


(21) 


can easily be computed. The sign of the additional terms 
is positive when path r used link (i,d) in the same direction 
as path p, and negative otherwise. With these components 
one can define a vector change in path flow, AQ such that 

1P('q + j3 Z?Q) ^ IP('q) (22) 


is a step size, which is selected by an appropriate one- 
dimensional search procedure. 

Several loadings should be considered in the 
design. Maximum hourly flows during the day and fire fighting 
demands are nomally used, but often the low demands, as may 
occur at night, have to be considered as well. For each 
loading, an initial flow distribution has to be specified 
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which satisfies eoatiauity at all nodes. A number of head 
constraints (7) are written for each loading, and the 
entire set is included in a single liP matrix. Once the LP 
has been solved, the gradient section modifies the flow 
flistribution for each of the loadings, using the results of 
the LP. 

Since the initial flows for each loading are quite 
arbitrary , there may not exist a set of segment diameters 
such that the bead line constraints (17) are satisfied for all 
loadings. Therefore, two new variables are added in each of 
these constraints - one with a positive and the other with 
a negative sign, and both required to be non-negative. These 
variables are assigned a large penalty coefficient in the 
objective function and can therefore be viewed as artificial 
variables. Each such variable may be viewed as a dummy 
valve in the path, able to take up the excess head for the 
specified flow in the path. If a dummy valve does appear 
in the final optimal solution, this means that a real valve 
will have to be installed at the location specified, and 
operated accordingly. When an actual valve does exist in 
the system it is represented in the LP by the head loss it 
provides under each loading. 

The iterative LPG procedure is terminated when 
any of several criteria is met (no significant improvement 
from one iteration to the next, specified number of flow 
change iterations exceeded , etc ,) . 
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When a pump is to be designed, the head it has to 
add for each of the loadings is the decision variable. These 
variables are introduced into the constraints (7) and (17) 
with the proper sign. An iterative procedure is used in 
dealing with the non-linear cost vs, head function for pumps. 
For reservoirs, the decision variable is the elevation at 
which it is to be located. This elevation appears in all 
constraints for paths ending at the reservoir. A linear 
cost vs. elevation relation is used in the objective function, 

2.5.2 Hon-linear Programming Models 

54 - 

An early effort was made to use a gradient-like 

technique in optimizing the design of a pipe network under 

one loading. The original work, which was not published, 

constituted the basis for lemieux’s thesis . Pipe diameters 

are changed, one pipe at a time, according to the derivative 

of objective function with respect to pipe diameters. After 

each change in a diameter the new network is solved, ; using 

the Newt on-Raph son method, and the last Jacobian of this 

solution appears in the computation of the derivatives. 

29 

Pitchai formulated a non-linear integer programming 

problem in seeking the optimal diameters of a network, and 

solved it by combining random search and examination of 

20 

adjacent design points. Jacoby used a gradient-approxi- 
mation method to seek minimum of a merit function which 
combined the objective function and penalties for violating 
head and continuity constraints for a single loading. Since 
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the search is conducted in the region which is hydraulically 
infeasible (it is an exterior point method), if the search 
terminates prematurely one may not have a feasible hydraulic 
design. A more detailed analysis of Jacoby’s paper may be 
found in Reference 35. 

Cembrowicz and Harrington^ have dealt with 
minimization of the capital cost of a pipe network designed 
to operate under one loading. Using graph theory they 
claim to decompose the problem such that the non-convex 
objective function is broken into subsets of convex functions, 
Bach function relates to a pipe or loop and is minimized 
separately, using a method of feasible directions. The 
number of optimizations may be very large and they have to 
be scanned to locate the global optimum, so that computa- 
tionally the method does not seem practical. 

Watanatada^ has developed a method for optimal 
design of pipeline networks supplied at a number of nodes, 
and applied it to real networks of moderate size. Bor a 


network with P pipes and M nodes, 

of 

which MS are supply 


nodes, the problem is 




P 


MS 

(23) 

Min Cj(D, H, Q) = 

+ 

-2 — St- , 
k=1 ^ 

subject to H, Q) = 0, 

k 

" 1 , * . . . , 

(24) 

Up ^ MIN, 

P 

= 1,,..,,P 

(25) 


k 

= 1 . , . ,M : 

(26) 
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- ^ 0, k = 1. ,MS (27) 

where D = a vector of P pipe diameters, 

H = a vector of M node heads, 

Q = a vector of MS supply rates at supply nodes. 

Up “ Up(Pp), the cost per unit length of pipe as 
a function of its diameter, 

“ ^k^^k’ ^k^* cost of supplying at a 
head 

QRjj. = algebraic sum of flows leaving the node. 

This constrained optimization formulation was converted into 
an unconstrained one by using a variable transformation due 
to Box^. Utility variables Z^, i = 1 (P + M + MS) 
are defined by 


Dp = DMIU + Zp , 

p “ 1,,,,*,P 

(28) 

% = . 

k = 1 , . . . . ,M 

(29) 

QNk = - , 

k = 1 , , . . ,MS 

(50) 

and the problem now becomes 



Min Gp(Z) 


(51) 

subject to QRj^(Z) = 0, 

k — 1 , , • o » ,M 

(52) 


A new function is now defined by combining const- 
raints (52) with fche objective function (51) according to 
the method suggested by Haarhoff and Buys ; 

M ^ o ' ''1 ' 

Minfp^CZ) = o’^(Z) * s 
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The subscript r is an iteration counter; and V/ are 
penalty multipliers - the first is updated at each iteration 
and the second is a preassigned fixed constant. The variable 
metric method of Fletcher and Powell^ ^ was used to minimize (53) j 
and was found to be superior to the Fletcher and Eeeves method 
of conjugate directions even though the former required 
more computer memory, Watanatada’s method has the advantage 
that the flow solution is incorporated directly into the 
optimization, and one therefore does not need a network 
solver as a separate program. The danger is, however, that 
if the procedure terminates promaturely the solution may not 
be feasible hydrolically. 

1 

A method based on Abadic’s G-RG method was 
developed by Shamir^ , in which the design and operation 
under a nimiber of loadings are to be optimized. For a network 
with N nodes , operating under 1 loadings the problem is s 

Min F(d, u, X, s) = f(d) + 21 w\^(d, u^, x^, s^) , (34) 

1=1 

subject to d ^ I) (35) 

u^ £ , V 1 (36) 

G^(d, u£ s^) = 0 \f 1 ^37} 

x^ = [”x/ G^(d, u^, x^, s^)*^ = 0 £ X (38) 

where d = the design variables (pipe diameters, pump capa- 
cities), which have to belong to the set D. 
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u = the operation variables (valves and pumps 01f/0F!P) 

for the 1 loading, which have to belong to the set u . 

X = the dependent variables (heads, consumptions) of 
th 

the 1 flow solution, which have to be within ranges given 

by 

= the independent (fixed) variable of the l"^^ flow 

solution. 

f = cost function of the design. 

= cost of operation for the 1^^ loading, 
w^ = weights, 

0 = a set of simultaneous node continuity equations 
(eq, (1)) for the 1^^ loading. The individual equations 
are = 0, j = 1 , . , . . ,N. 

It is assumed that the sets b, and X^ simply 
specify a range of values for the corresponding variables. 
The L sets of equations (37) are combined with the objective 
function (34) to form the Lagrangians 


1 H 

c£(d, u, X, s, A) = ^(d, u, X, s) + Z. ^ 

1=1 j=1 


1 „ 1 ,^ 1 1 
^G^(d, u , X , 

ti 


= S' • W ^59) 

where T stands for the transpose. At any point (d, u) the 
following has to hold 

Since G and x can be separated into independent 
flow problems ♦ eqn. (40) can be decomposed into 
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1 = 1 , . . . ,L (41 ) 

The Langrange multipliers are therefore solved in L groups 
of N values each from eqn. (41). The matrix psV^x^jis 
the last Jacobian of the 1^^ flow solution by the Newton-' 
Eaphson method » and is therefore available directly as a 
by-product of the flow solution. Once the 's have been 
obtained, the components of the (reduced) gradient are 



At a (local) optimum F = 0. At any other point - F 

points in the direction of the steepest descent of F, while 

changes in flow solutions due to a move in this direction 

are already taken into account. A move in this direction is 

now made, using a one-dimensional search procedure. At the 

new point the L flow problems are solved and the reduced 

1 

gradient calculation repeated. The constraints on d, u , 
and x^ are used in selecting the step size in each move. 

The search is terminated by given criteria for improvement 
between iterations, value of the g 2 ?adienfc , number of 
iterations etc. 
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Because flow solutions are computed at each step, 
and because the constraints on the decision variables are 
not violated during the moves, the current solution is 
feasible (this is essentially an interior point method), 
and if it terminates prematurely, one at least has a feasible 
solution - -^stiich is better than the one having started with. 

2.5.3 Other Methods 

7 

Deb and Sarkar based their method on the concept 
of equivalent diameters, for least cost design of a network 
operating under a single leading and in which the heads were 
assumed to be known. The results are therefore quite limited 

in application. 

6 

Deb again used a similar combination of the flow 
and cost equation, and developed a method for optimal design 
of a system consisting of a pumping station, an elevated 
reservoir and a pipe network fed from it by gravity. The 
results are further restricted by the fact that the shape 
s£ pressure surface over the network is assumed to have a 
specific form. 

2 .6 Operation Over Time 

Some of the methods described above can be used 
to reach optimal operating decisions for existing networks. 
This is done by fixing all design variables at their actual 
values and carrying out the optimization for the operational 
variables. Special methods have, therefore, to be developed 
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for optimization of operation over time, in wiiich the only 

decision variables are the operation of pumps and the 

setting of valves during the specified time horizon, 

1 1 

Dreizin et al. used a hydraulic simulator of a 
particular water system as the basic building block in a 
program which attempted to improve operating policies. The 
decision variables were those water levels (called set- 
points) in specified reservoirs at which pumps are to be 
switched ON or OT'i'. No algorithmic optimization was found 

i ■ 

'for solving the problem, and a sequence of simulations with 
response surface analysis (a gradient-like search) was used. 

o 

Some work in the City of Philadelphia resulted 
in selection of operating policies over a day based on a 
comparison of costs for several proposed policies. No 
optimization was attempted. 

Sterling and Coulbeck^® * optimized pumping costs 
in a watei* system using dynamic programming and a two level 
hierarchical approach. A group at the University of 
Cambridge, England, has been engaged in development of on- 
line control and optimization of the operation of regional 
water systems^ which are being implemented and tested 
in the field , 
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3. MODEL FOMTEATION 


The model presented here is developed "by Alperovits 

2 

and Shamir . The model uses Linear Progi^mming Gradient 
(LPG) method to optimize the solution and is the first to 
incorporate flow solution into the optimization procedure, 
without making any assumptions about the nature of pressure 
surface as attempted by the earlier investigators. 

The Linear Programming Gradient method deals with 
loop^ networks and decomposes the optimization problem 
into a hierarchy of 2 levels as depicted in Figure 1 . 



Figure 1 . Schematic Diagram of the LPG Method, 
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The first step in developing the LP& method is 
to consider optimization of the design when the distribution 
of flows in the network is assumed to be known. We adopt 
the formulation given by objective function (8) and constraints (5 
(7)*(J7)in which the lengths of the segments of candidate 
diameter in each link are the decision variables, For any 
Q the optimal cost of the network may be written as 


Cost = LP ( Q ) 

The next stage is to develop a method for syste- 
matically changing Q with the aim of improving cost. The 
method for changing Q is based on the use of the dual vari- 
ables, which aid in defining a gradient move. As the flow 
distribution should always be balanced at the nodes hence 
equal corrections (sign of the correction depending on the 

flow direction) are applied to each pipe in a given loop. 

2 

Shamir proposed that the component of the gradient vector 
G corresponding to loop p should be as 


(Cost) 




(43) 


Summations are performed only for the appropriate links, 
i.e,, those belonging to the p loop. 

Quindry et .al . modified the gradient vector G 

taking into consideration the interactions of paths with 
each other as well and proposed their formulation for G^ as 
given in equation (19). Their formulation was found to give 



better results compared to equation (45). The proof for 
the equation (19) is given in Appendix 1. 

^ IflcoJ^Porating Pumps in the Formulation 
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When there are to be ptmps in the system, the 
locations at which pumps may be installed are selected by 
the designer, but since program can set certain pump capa- 
cities to zero, if that is the optimal solution, the program 
actually selects the locations at which pumps will be 
installed . The decision variable associated with each 
location is the head to be added by the pump. If one 
denotes by XP(t) the head added by pump number t, then the 
head constraints of the type (7) for paths with pirnips become 

HMIN < H . XP(t) 1 SI £ J. .X ^ HMA2 (44) 

where the first summation is over the pumps in the path. 

For any path which has pumps in it, be it a closed loop or 
an open path, equation (1?) had to be modified in a similar 
manner, and it then becomes 


^ ^ ^ijm ^ijm - ^ 


(45) 


i ,3 m t 

In equations (44) and (45) the signs in front of the various 
terms depend on the direction of flow. 

The decision variables for the pumps XP(t) have 
to be introduced linearly into the objective function if 
the problem is to remain a linear program. This is done by 
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considering the cost of the ptanp as a linear functions of the 
capacity, i.e,, its rated horsepower. However, the cost 
per horse-power decreases with increasing pump capacity. 
Successive approximations are used in the program to cope 
with the non-linearity of the cost curve. The power needed 
to operate the pump is given hy 

HP = Y , Q , XP/y^ (46) 

where Y is a coefficient, Q is the flow, XP is the head 
added hy the pump, and is the efficiency. Por a fixed 
discharge and a fixed efficiency, equation (46) hecomes 

HP = K . XP (47) 

where K is a constant. 

The operating cost of the pump is therefore 
linearly proportional to the decision variable, which is XP. 
An iterative procedure is developed to take into account 
the non-linearity of the capital cost as follows? 

( 1 ) Assume values for the cost per HP for each pump 
location. 

( 2 ) Solve the linear program with these values as the 
coefficient of the XP in the objective function. 

( 3 ) For the resulting XP after the IP has been solved, 
compute the cost per HP . If all values are close to 
those assumed, this step is complete and one proceeds 
to a flow iteration by the gradient method. Otherwise 
one takes the new costs and solves the program again. 

This procedure is found to work well. ^ 
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^ Reservoirs xn the Forniulation. 

The decision variable for a reservoir is the 
elevation at which it is to be located. An initial elevation 
is assumed, then SR is the additional elevation where the 
reservoir is to be located, relative to its initially assumed 
elevation (2R. is addition to reservoir height and XR” is 
deduction from its height). Path equations have to be formed 
between the reservoir at node s and nodes in the network. 

For node n 

4 HO3 * XR; - JE; 1 ^ ^ X. / HMAX^ 

X , j m 

(48) 

where HOg is the initial elevation of the reservoir at node 

s, and ZRg is the additional elevation to be selected by 

the program. The coefficient of ZR„ in the objective function 

s 

is the cost of altering the location of the reservoir by 
1 unit (1m). 

3.3 Selection of Candidate Diameters 

At the outset, the list of candidate diameters 
for each link is based on a suitable hydraulic gradient say 
0.001. Out of the available pipe diameters, the diameter 
is selected which gives the nearest hydraulic gradient to 
the above value. Candidate list consists of 3 pipe 
diameters, one larger and one smaller than the one selected 


above . 
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If the first selection of candidate diameters 
results in an optimal LP solution, the list may have to he 
modified after the flows in the network have been changed 
by the gradient move. The modifications in the list of the 
candidate diameters are based on the following rules; 

(a) If in the optimal IP solution a link is made entirely 
of one diameter; then for the next IP, the list is made of 
5 diameters, the existing one and both its neighbours. 

(b) When in the optimal solution a link is made of 2 dia- 
meters, the list for the next LP is made of these two, plus 
one adjacent to that diameter which has the larger of the 
two lengths . 

Both cases result in a list of only 3 diameters 
for each link. In going from one LP solution to the other, 
the list for a particular element may remain unchanged, 
or a diameter may be dropped off one end of the list and a 
new one added at the other. 

3.4 Overview of the Algorithm 

(a) The first step is to decide about the initial flow 
distribution knowing the demand at each of the nodes. 

(b) Objective function and constraints are foimulated as 
described previously. 

(c) The LP is solved to give the optimum results for the 
initial flow distribution. The gradient vector is 
calculated knowing the dual variables of the LP. 


Cd) The magnitude of the gradient vector is calculated as 
follows 

MASHITODE = (af + 8^ * (49), 

where G-^ , *j2’****’^n components of gradient vector. 
If the magnitude is less than a prespecified value the 
problem is solved. If not so, the flows are to be 
changed in the gmdient direction. 

(e) The gradient vector is normalized by dividing every 
component by MAGNITUDE, 

(f; The gradient vector is multiplied by a trial step 

length to give the required flow corrections in loops. 

(g) After applying these corrections, if any of the flow 
direction changes, the step length is reduced by a 
factor (say half) and new flow corrections are 
calculated . Program is terminated if the step length 
becomes less than a prespecified value. 

(h) Eor the new flow pattern, the LP is formulated and 
solved. 

(i) If the optimum cost of the new IP comes out to be 
less than the preceding IP, the new gradient vector 
is calculated and control is shifted to step (d). 

If the optimum of the new IP is more than the preceding one 
then the step length is reduced by a factor and the preceding 
gradient move and flow pattern is used to get the new flow 
pattern and control is shifted to step (h). If the step 
length becomes less than a prespecified value, the program 
is teimiinated and the present solution is taken to be the 
optimum,: 



40 


Qg WATER DISTRIBUTION SYSTEM USING THE fROGHAM 

A computer program in Fortran has 'been written 
based on the LPG method as given in the last section. The 
program can be used to design a new water distribution system 
given the consiimption pattern during peai hours at various 
nodes of the system. Some salient features regarding the 
use of the program are as given below s 

i) The critical consumption pattern, i.e, the peaic 
demand should be known. 

ii) An initial flow distribution along the pipes is to 
be given as input data. The given flow distribution 
should satisfy the continuity of flow requirement at 
all nodes of the system. 

iii) In case of more than one source of supply, the amount 
of water to be withdrawn from different sources 
should be decided beforehand, i.e. the program does 
not take into consideration the optimal allocation 
of resources. 

iv) The suitable positions for the booster pumps should 
be given as input data. If a pump is required at a 
specified position to minimize the cost, the 
program will give the output in terms of the head to 
b^^ at that point. 
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Example Problems 


The program is tested for a set of problems with 
increasing network size. Execution time for different 
problems is reported for Dec- 10 system (available at IIT/ 

Kanpur), Details of problems along with solutions are as 
follows t 


The set of available diameters and the cost of 
pipes and pumps are common to all problems. The cost of 
pumps is taken to be Es. 1000/- per installed horse-power. 
The cost data for pipes are given in Table 1 . The cost for 
the source of supply is to be given in terms of the increase 
in cost for 1 m increase in head at the supply node. An 
initial head is assumed at the supply node and the program 
gives the increase or decrease in head required to minimize 
the cost. As the lives of various components in the distri- 
bution system are different so all the capital costs are 
converted to annual costs and the total annual cost of the 
system is minimized. 

Capital costs can be converted to annual costs by 
multiplying by the Capital Hecovery Eactor (OSE), 


CEE 


- id i) 


n 


(1 


+ i)^ - 1 


in which i is the rate oi and n is the estimated 

life of the component . Assuming 10^ rate of interest and 
15 years and 40 years to be the lives of pumps and pipes 
uespeotively the factors come out to be; 
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CRF (Pump) = 

(1 + 0 . 1 )'^ - 1 


0.13147 


CHP (Pipes) = 

♦ 


0.10225 


To calculate the increase in cost due to increase 

of supply head of pump hy 1 m, the following procedure is 
adopted s 


Assuming the pump cost is Rs. 1000/- per installed 
horse -powers 

Hence the increase in cost due to 1 m increase in head 
= 1000 X 2.7323 X 10“"^ X Q = 0.27323 X Q Rs/m 


where Q is water supplied by pump in 1pm. 

Assuming that the pump operates at this rate for 8 hours a 
day, the increase in annual cost due to operation 

= Q ^ Q • . 23 ^ 0.14873 Q Rs/m 


where cost of power = Rs, 0.25/lcwh and efficiency of pump = 0.8. 
Hence increase in cost due to 1 m increase in supply head 
= (Capital cost x CEB’) + Operation cost 
= [( 0.27323 X 0 . 13147 ) + 0.1487^ Q 
= 0.18465 xQ Rs/m 

Problems 1 to 8 r Distribution networks for the problems are 
given in the following figures; 

Problem 1 - Figure 2 
Problem 2 - Figure 3 
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Problem 4 

Figure 5 

Problem 5 

Figure 6 

Problem 6 

Figure 7 

Problem 7 

Figure 8 

Problem 8 

Figure 9 

Every problem has a pump serving as the source and 

a booster pump, Node no. 1 is supply node in each problem 

and rest all nodes 

have 100 1pm demand. Hence supply at 

node 1 is given as 

follows : 

Problem 1 

600 1 pm 

Problem 2 

800 1 pm 

Problem 3 

- 900 1pm 

Problem 4 

“ 1100 1 pm 

Problem 5 

- 1 200 1 pm 

Problem 6 

- 1 400 1 pm 

Problem 7 

- 1500 1pm 

Problem 8 

- 1600 1 pm 

The minimum pressure allowed at each node is taten to be 

15 m. 


The location of booster pump to be provided in 

each problem is given as follows; 

Problem 1 

Pipe No. 2 

Problem 2 

pipe No, 6 

Problem 5 

- Pipe No. 6 

Problem 4 

- Pipe No. 6 

Problem 5 

- Pipe No, 6 

Problem 6 

Pipe No. 6 

Problem 7 

- Pipe No. 6 

problem 8 

- Pipe No. 9 
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The initial flow distribution and candidate 


diameters for each 

pipe 

are given in the following tables 

Problem 1 

- 

Table 2 

Problem 2 

- 

Table 3 

Problem 3 

— 

Table 4 

Problem 4 

- 

Table 5 

Problem 5 

- 

Table 6 

Problem 6 

- 

Table 7 

Problem 7 

- 

Table 8 

Problem 8 

- 

Table 9 

The structures 

of initial linear programmes are 

given in following 

tables : 

Problem 1 

- 

Table 10 

Problem 2 

- 

Table 1 1 

Problem 3 

- 

Table 12 

Problem 4 

- 

Table 13 

Problem 5 


Table 14 

Problem 6 

- 

Table 15 

Problem 7 


Table 16 

Problem 8 

— 

Table 17 


Increase in cost due to 1 m increase in supply- 
head for various problems is as following; 


Problem 1 
Problem 2 
Problem 3 
Problem 4 


1 10 .79 Rs/m 
147.67 Rs/m 
166.18 Es/m 
203.11 Es/m 



Problem 5 
Problem 6 
Problem 7 
Problem 8 


221.58 Rs/m 
258.51 Rs/m 
- 276.97 Rs/m 

~ 295.44 Rs/m 
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Table 1 

Available P ipes , Their Costs, and Hazen-Williaa Coefficients 

CHI' for pipes = 0.10225 


S.No. 

t 

I Pipe 

J diameter 

I (mm) 

1 

Cost 

Rs/m 

— 5 

I 

I Annual cost 

5 Rs/m 

f 

t 

f 

1 

1 

1 Hazen-Willian 

I coefficient 

f 

1 

I 

1 

15.0 

2.00 

0.2045 

» 

140.0 

2 

20.0 

3.00 

0.30675 

140.0 

3 

25.0 

4.50 

0.46012 

140.0 

4 

32.0 

5.00 

0.51125 

140.0 

5 

40.0 

7.00 

0.71575 

140.0 

6 

50.0 

8.00 

0.818 

140.0 

7 

63.0 

1 1 .00 

1.12475 

140.0 

8 

80.0 

15.50 

1 .58487 

140.0 

9 

100.0 

22.00 

2.2495 

140.0 

10 

125.0 

, 35.50 

3.62987 

140.0 

11 

150.0 

47.00 

4.80575 

140.0 

12 

200.0 

83.00 

8.48675 

150.0 

13 

250.0 

112,00 

11.452 

130.0 

14 

300.0 

145.00 

14.82625 

130.0 

15 

350.0 

187.00 

19.12075 

130.0 

16 

400.0 

221 .00 

22.59725 

130.0 

17 

450.0 

278.00 

28,4255 

130.0 

18 

500.0 

329,00 

33.64025 

130.0 

19 

600.0 

463.00 

44.34175 

130.0 


Note ; 

1) According to the present practice in U.P. Jal Nigam, upto 
150 caa diameter, P?C pipes are provided and beyond this 
size, AC pipes are used. The costs and Hazen-WilliaB 
coefficients are taken accordingly. 

2 ) 15 pip® bas also been taken in the set of available 
pipes in the present case but usually such small pipes; 
should not be considered "while designii^ the distribution 
system. The optimum results will come out to be different 
if saall pipes are not included in the list of available 
pipes. 


Table 2 

Pipe Data for Problem - 1 


Pipe 

No. 

length 

(a) 

r 

Initial Plow 
(1pm) 



f 

Candidate diameters 

I (mm) 

f 

t 

1 

1000.0 

600*0 

150, 

200, 250 

2 

1000.0 

280.0 

100, 

125, 150 

3 

1000,0 

180.0 

80, 

100, 125 

4 

1000.0 

10.0 

32, 

40, 50 

5 

1000.0 

220.0 

100, 

125. 150 

6 

1000.0 

90.0 

63. 

80,^100 

7 

1000.0 

10.0 

32, 

40, 50 

8 

1000.0 

110.0 

80, 

100, 125 


Table 3 

Pipe Data for Problem - 2 

I 

Pipe ; 
No. ; 

t 

f 

Length 

(a) 

1 — 

1 Initial Plow 
; (1pm) 

f 

I 

f 

1 Candidate diameters 

I (mm) 

f 

f 

1 

1000.0 

800.0 

150, 

200, 

250 

2 

1000.0 

590.0 

150, 

200, 

250 

3 

1000.0 

210.0 

100, 

125, 

150 

4 

1000.0 

10.0 

32, 

40, 

50 

5 

1000.0 

110.0 

80, 

100, 

125 

6 

1000-0 

280.0 

100, 

125, 

150 

7 

1000.0 

180.0 

80, 

100, 

125 

8 

1000.0 

10.0 

32, 

40, 

50 

9 

1000.0 

90.0 

63, 

80, 

100 

10 

1000.0 

10.0 

32, 

40, 

50 

11 

1000.0 

110.0 

80, 

100, 

125 



TaTale 4 

Pipe Data for Problem - 3 


k 

Pipe ; 
No. : 

f 

1 

Length 

(m) 

1 ' ' ■ 

; Initial flow 

; (Ipm) 

1 

f 

' ' ' " 

t 

; Candidate diameters 

1 (mm) 

f 

f 

1 

1000.0 

900.0 

150, 

200, 

250 

2 

1000.0 

580.0 

150, 

200, 

250 

3 

1000.0 

200.0 

100, 

125, 

150 

4 

1000.0 

10.0 

32, 

40, 

50 

5 

1000.0 

220.0 

100, 

125, 

150 

6 

1000.0 

280.0 

100, 

125, 

150 

7 

1000.0 

180.0 

100, 

125, 

150 

8 

1000.0 

10.0 

32, 

40, 

50 

9 

1000.0 

90.0 

63, 

80, 

100 

10 

1000.0 

10.0 

32, 

40, 

50 

11 

1000.0 

100.0 

80, 

100, 

125 

12 

1000.0 

10.0 

32, 

40, 

50 

13 

1000.0 

110.0 

80, 

100, 

125 
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Table 5 

Pipe Data for Problem - 4 


Pipe 

Ho. 


Initial flow 
(1pm) 

“T — 

t 

; Candidate diameters 
; (mm) 

t 

f 

1 

1000.0 

1100.0 

200, 

250, 

300 

2 

1000.0 

780.0 

150, 

200, 

250 

3 

1000.0 

200.0 

100, 

125, 

150 

4 

1000.0 

10.0 

32, 

40, 

50 

5 

1000.0 

220.0 

TOO, 

125, 

150 

6 

1000.0 

480.0 

125, 

150, 

200 

7 

1000.0 

190.0 

100, 

125, 

150 

8 

1000.0 

10.0 

32, 

40, 

50 

9 

1000.0 

190.0 

100, 

125, 

150 

10 

1000.0 

90.0 

63, 

80, 

100 

11 

1000.0 

10.0 

32, 

40, 

50 

12 

1000.0 

90.0 

63, 

80, 

100 

13 

1000.0 

10.0 

32, 

40, 

50 

14 

1000.0 

100.0 

63, 

80, 

100 

15 

1000.0 

10.0 

32, 

40, 

50 

16 

1000.0 

110,0 

80, 

100, 

1 25 




Table 6 


Pipe Data for Problaa - 5 


Pipe ; 
No. ; 

Length 

(m) 

— J 

t 

Initial flow 

I (1pm) 

t 



J Candidate diameters 

! Cnmi) 

1 

1 

1 

1000.0 

1200.0 

200, 

250, 

300 

2 

1000.0 

880.0 

150, 

200, 

250 

5 

1000.0 

200.0 

100, 

125, 

150 

4 

1000.0 

10.0 

32, 

40, 

50 

5 

1000.0 

220,0 

100, 

125, 

150 

6 

1000.0 

580.0 

150, 

200, 

250 

7 

1000.0 

200.0 

100, 

125, 

150 

8 

1000.0 

10.0 

32, 

40, 

50 

9 

1000.0 

280.0 

100, 

125, 

150 

10 

1000.0 

180.0 

80, 

o 

o 

125 

11 

1000.0 

10.0 

32, 

40, 

50 

12 

1000.0 

90.0 

63, 

80, 

100 

15 

1000.0 

10.0 

32, 

40, 

50 

14 

1000.0 

100.0 

63 5 

80, 

100 

45 

1000.0 

10.0 

32, 

40, 

50 

16 

1000.0 

100.0 

63, 

80, 

100 

17 

1000.0 

10.0 

32 , 

40, 

50 

18 

1000.0 

110.0 

80, 

100, 

125 
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Table 7 

Pipe Data for Problem - 6 


' ' ' " f 

t 

Pipe \ 

No. ; 

1 

f 

length 

(m) 

t 

t 

! Initial flow 

1 (1pm) 

! 

f 

i 

f 

Candidate diameters 
! (mm) 

t 

1 

1000.0 

1400.0 

200, 250, 300 

2 

1000.0 

1080.0 

150, 200, 250 

3 

1000.0 

200.0 

100, 125, 150 

4 

1000.0 

10.0 

32, 40, 50 

5 

1000.0 

220.0 

ioo , 125, 150 

6 

1000-0 

780.0 

150, 200, 250 

7 

1000.0 

310.0 

125, 150, 200 

8 

1000,0 

10.0 

32, 40, 50 

9 

1000.0 

370.0 

125, 150, 200 

10 

1000.0 

270.0 

100, 125, 150 

11 

1000.0 

10.0 

32, 40, 50 

12 

1000.0 

180.0 

80, 100, 125 

13 

1000.0 

10.0 

32, 40, 50 

14 

1000.0 

210.0 

100, 125, 150 

15 

1000.0 

10.0 

32, 40, 50 

16 

1000.0 

100.0 

80, 100, 125 

17 

1000.0 

10.0 

32, 40, 50 

18 

1000.0 

110.0 

80, 100, 125 

19 

1000.0 

90.0 

63, 80, 100 

20 

1000.0 

10.0 

32, 40, 50 

21 

1000.0 

110.0 

80, 100, 125 


Table 8 

Pipe j)ata for Problem - 7 


T 

Pipe i 
No. ; 

f 

Length 

(m) 

-7 — 

f 

1 Initial flow 
j (1pm) 

t 

"T" — ■ 

t 

J Candidate diameters 

1 (nrni) 

t 

t 

1 

1000.0 

1500.0 

200, 250, 300 

2 

1000.0 

1180.0 

200, 250, 300 

3 

1000.0 

310.0 

100, 125, 150 

4 

1000.0 

10.0 

32, 40, 50 

5 

1000. 0 

220.0 

100, 125, 150 

6 

1000.0 

770.0 

150, 200, 250 

7 

1000.0 

300.0 

100, 125, 150 

8 

1000.0 

10.0 

32 , 40 , 50 

9 

1000.0 

370.0 

125, 150, 200 

10 

1000.0 

270.0 

100, 125, 150 

11 

1000.0 

10.0 

32, 40, 50 

12 

1000.0 

180.0 

80, 100, 125 

13 

1000.0 

10.0 

32, 40, 50 

14 

1000.0 

200.0 

100, 125, 150 

15 

1000.0 

10.0 

32, 40, 50 

16 

1000.0 

210.0 

100, 125, 150 

17 

tooo.o 

10.0 

32, 40, 50 

18 

1000.0 

110.0 

80, 100, 125 

19 

1000.0 

90.0 

63, 80, 100 

20 

1000.0 

10.0 

32, 40, 50 

21 

1000.0 

100.0 

.63, 80,' 100 

22 

1000.0 

10.0 

32, 40, 50 

23 

1000.0 

110.0 

80, TOO, 125 
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Table 9 

Pipe Data for Problem - 8 


I 

1 

Pipe ; 

Length 

t ■ ■ 

f 

I Initial flow 

t ' " ' ' ' " 

* 

; Candidate diameters 

No. I 

t 

f 

(m) 

' (1pm) 

t 

1 

f 

I 

f 

(mm) 


1 

1000.0 

1600.0 

200, 

250, 

300 

2 

1000.0 

1170.0 

200, 

250, 

300 

5 

1000.0 

300.0 

125, 

150, 

200 

4- 

1000.0 

10.0 

32, 

40, 

50 

5 

1000.0 

330.0 

125, 

150, 

200 

6 

. 1000.0 

770.0 

150, 

200, 

250 

7 

1000.0 

300.0 

125, 

150, 

200 

8 

1000.0 

10.0 

32, 

40, 

50 

9 

1000.0 

370.0 

125, 

150, 

200 

10 

1000 .0 

270.0 

100, 

125, 

150 

11 

1000.0 

10.0 

32, 

40, 

50 

12 

1000.0 

180.0 

80, 

100, 

125 

13 

1000.0 

10.0 

52, 

40, 

50 

U 

1000.0 

200.0 

100, 

125, 

150 

15 

1000.0 

10.0 

32, 

40, 

50 

16 

1000.0 

200.0 

100, 

125, 

150 

17 

1000.0 

10.0 

32, 

40, 

50 

18 

1000.0 

220.0 

100, 

125, 

150 

19 

1000.0 

90.0 

63, 

80, 

100 

20 

1000.0 

10.0 

32, 

40, 

50 

21 

1000.0 

100.0 

80, 

100, 

125 

22 

1000.0 

10.0 

•y 

,, 5 

40, 

50 

23 

1000.0 

100.0 

80, 

100, 

125 

24 

1 000 .0 

10.0 

32, 

40, 

50 

25 

1000.0 

110.0 

80, 

too. 

125 
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Table 10 

Initial linear Program for Problem - 1 


t 

Begin i 
Node ; 

f 

t 

End 

Node 

— r — ' — — — 

f, 

; Number of sections connected between 

\ the nodes 

1 

t 



Presstrre Eauations 

1 

6 

5, 8 

1 

7 

U 2, 3, 6 



Loop Eauations 

2 

2 

2, 3. -4, -5 

4 

4 

4, 6, -7, -8 


Table 11 

Initial Linear Program for Problem - 2 

f 

Begin I 

End 

t 

1 

; Numbers of sections connected between 

Node ; 

• 

, f 

Node 

; the nodes 



Pressure Eauations 

1 

5 

1 , 5 

1 

8 

1, 2, 3, 11 

1 

9 

1, 2, 6, 7, 9 



Loop Equations 

2 

2 

2, 3, -4, -5 

3 

3 

6, 7, -8, -3 

6 

6 

8, 9, -10, -11 
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Table 12 

Initial linear Program for Problem - 3 


% 

Begin \ 
Node ; 

1 

End 

Node 

— r — — 

1 ■ 

1 Numbers of pipes connected between 

I the two nodes 

t 



Pressure Eauations 

1 

8 

U 5, 13 

1 

9 

1t 2, 3» 11 

1 

TO 

U 2, 6, 7, 9 



Loop Equations 

2 

2 

2, 3, -4, -5 

3 

3 

6, 7, -8, -3 

6 

6 

8, 9, -10, -11 

5 

5 

4, 11, -12, -13 


Table 13 

Initial Linear Program for Problem -4 

^ -P 

f 

Begin 1 

End 

t 

; Numbers of pipes connected between 

Node : 

■ ■ !'■ 

f 

Node 

I the two nodes 

j 

f 



Pressure Equations 

1 

10 

1, 5, 16 

1 

11 

1, 2, 3, 14 

1 

1 2 

1,2,6,7,12 

1 

9 

1,2,6,9,10 



Loop Equations 

2 

2 

2, 3, -4, -5 

3 

3 

6, 7, -8, -3 

4 

4 

9, 10, -11, -7 

7 

7 

8, 12,-13, -14 

6 

6 

4, 14, -15, -16 




Table 14 

Initial Linear Program for Problem - 5 




T — 

t 

Begin I 

Node 1 

1 

l 

End 

Node 

— r — - 

! Numbers of pipes connected between 

I the two nodes 

f 



Pressure Equations 

1 ■ 

10 

1, 5, 18 

1 

11 

1,2,3, 16 

1 

12 

1f256f7f^4 

1 

13 

1, 2, 6, 9, 10, 12 



loop Equations 

2 

2 

2, 3., -4, -5 

3 

3 

6, 7, ”8, -3 

4 

4 

9, 10, -11, -7 

8 

8 

11, 12, -13, “14 

7 

7 

8, 14, -15, -16 

6 

6 

4, 16, -17, “18 



Initial 

Table 15 

Linear Program for Problem - 6 

T — 

t 

Begin I 
Node : 

f 

t 

End 

Node 

t 

1 Numbers of pipes connected between 
\ the two nodes 

i 

f ' ■ : - - 



Pressure Equations 

1 

10 

1,5, 18 

1 

11 

1,2, 3, 16 

1 

14 

1 , 2, 6, 7, 14, 21 

1 

15 

1, 2, 6, 9, 10, 12, 19 



Loop Equations 

2 

2 

2, 3, “4, “5 

5 

3 

6, 7, -8,-3 

4 

4 

9, 10, -11, -7 

8 

8 

11, 12,-13, -14 

7 

7 

8, 14, “15. “16 

6 

6 

4, 16, -17, “18 

12 

12 

13, 19, “20, -21 
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Table 16 

Initial linear Program for Problem - 7 


Begin | 

Node ; 

■ 1 

Snd 

Node 

,* Numbers of pipes connected between 

I the two nodes 

1 

10 

Presstare Bauations 

1, 5, 18 

1 

14 

1, 2. 3, 16. 23 

1 

15 

1 , 2, 6, 7, 14, 21 

1 

16 

1,2, 6, 9, 10, 12, 19 

2 

2 

Loop Eauations 

2, 3, “4, -5 

3 

3 

6, 7, -8, -3 

4 

4 

9, 10, -11, -7 

8 

8 

11 , 12, -13, -14 

7 

7 

8, 14, 15, -16 

6 

6 

4, 16, -17, -18 

12 

12 

13, 19, -20, -21 

11 

1 1 

15,21.-22,-23 



Initial 

Table 17 

Linear Program for Problem - 8 

■' ' . I 

Begin ; 

End 

, rT-,-,,-.- ri, - 

i Numbers of pipes connected between 

t 

Node 

; the two nodes • 

I 



Pressure Eauations 

; T ■ : ’ 

14 

1, 5, 18, 25 

1 

15 

1,2, 5, 16, 23 

1 

16 

1, 2, 6, 7, 14, 21 

; : : 1 ; 

17 

1, 2, 6, 9, 10, 12, 19 



L^^^ Equations 

2 


2, :3, -4, -5 

3 

3 

6, 7, -8, -3 

4 

4 

9, 10, -11, -7 

6 

6 

4, 16, -17, -18 

7 

7 

8. 14, -15, -16 

8 

8 

11, 12, -13, -14 

12 

12 

13, 19,-20,-21 

11 

11 

15, 21, -22, -23 

10 

10 

17, 23,-24, -25 
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5. BES'DLTS AED SISCUSSIQW 


The optimum designs obtained from execution of 
the program are given in following tables. 


Problem 1 

Table 

18 

Problem 2 

Table 

19 

Problem 3 

Table 

20 

Problem 4 

Table 

21 

Problem 5 

Table 

22 

Problem 6 

Table 

23 

Problem 7 

Table 

24 

Problem 8 

Table 

25 


Table 18 


Optimum Design for Problem - 1 


Pipe 

Plow 

Length 

(mi 

I 

; Segment 1 

f 

i ' ' ^ ' 

; Segment 2 

f 

No. 

(1pm) 

f 

roianeter 

1 (nm) 

t 

Length 

(nj 

f 

^Diameter 

1 (nm) 

f 

Length 

(ml 

1 

600.00 

1000.0 

100 

1000.00 

mm 

. 

2 

285.53 

1000.0 

80 

1000.00 

mm 

- 

3 

183.53 

1000.0 

80 

1000.00 

mm 

■ 

4 

2.57 

1000.0 

15 

1000.00 

- 

-- 

5 

216.47 

1000.0 

80 

799.05 

63 

200*95 

6 

86.10 

1000.0 

63 

1000.00 

- 

- , 

7 

13.90 

1000.0 

32 

1000.00 



8 

113.90 

1000.0 

65 

994.09 

50 

■5.91:,; .'i 


Optimum pumping head at Node 1 =35 + 18,26 = 53.26 m 
Minimm cost of the system = 11898.25 Rs /year 
Minimum cost for the initial flow pattern = 12095.24 Rs/year 
At pipe no. 2, no booster pump is req^uired. 
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Table 19 

Optimum Design for Problem - 2 


Pipe 

Wo . 

r 

Plow 

(1pm) 

"1 

t 

1 Length 

1 (m) 

t 

f 

1 

t 

1 Segment 

r * 

1 

f 

Segment 2 

i ' i 

J Diameter \ 

; (mm) ; 

f t 

— ^ 1 

Length Hliameter 
(mj ; (mm) 

f 

Length 

(m) 

1 

800.00 

1000.0 

125 

1000.00 

- 


2 

599.99 

1000.0 

125 

565.94 

100 

434.06 

5 

211 .77 

1000.0 

80 

996 . 64 

63 

3.36 

4 

0.01 

1000.0 

15 

1000.00 

- 

- 

5 

100.00 

1000.0 

50 

1000.00 

- 

- 

6 

288.22 

1000.0 

too 

1000.00 

- 

... 

7 

188.22 

1000.0 

80 

1000.00 

- 


8 

11 .50 

■ 1000.0 

52 

1000.00 

- 

- 

9 

99.72 

1000.0 

80 

667.62 

65 

322.38 

10 

0.27 

1000.0 

15 

1000.00 

- 


11 

100.27 

1000.0 

65 

1000.00 

mm 



Optimum pumping head at node 1 = 55 + '•2,6 = 47.6 m 

Minimum cost of the system = 18258.60 Hs /year 

Minimum cost for the initial flow pattern = 18842.00 Rs /year 

k% pipe no. 6, no booster pump is required. 
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Table 20 

Optidal Design for Problem - 3 


T ! r 


Pipe 

Plow 

; Length 

; Segment 1 

; Segment 2 

f ■ 

Ho. ; 

f 

1 

1 

f 

(1pm) 

J Cm) 

• i 

3)iameter I Len^h 
; (mm) 1 (m) 

f t 

1 

JDiameter 

; (mm) 

1 

Len^h 

(m) 

1 

900.00 

1000.00 

125 

1000.00 


- 

2 

' 584.18 

1000.00 

125 

304.82 

100 

695.18 

3 

199.94 

1000.00 

80 

1000.00 

- 

- 

4 

2.18 

1000.00 

15 

1000.00 

- 

- 

5 

215.81 

1000.00 

80 

1000.00 

- 

- 

6 

280.06 

1000.00 

100 

992.55 

80 

7.45 

7 

180.06 

1000.00 

80 

1000.00 



8 

10.73 

1000.00 

32 

913.37 

25 

86*65 

9 

90.79 

1000.00 

63 

1000.00 


- 

10 

9.20 

1000.00 

32 

1000.00 

- 

- 

11 

95.58 

1000.00 

63 

941 .89 

50 

58*11 

12 

13.63 

1000.00 

52 

1000.00 


- 

13 

113.63 

1000.00 

63 

1000.00 

- 



Optimum pumping head at node 1 = 35 18,22 = 53.22 m 

Minimum cost of the system = 2141T.4'3 Rs/year 

Minimum cost for the initial flow distribution = 21648.22 Rs/year 
Ho booster pump is required at pipe Ho. 6. 
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Table 22 

Optiiam Design for Problem -- 5 


T 


Pipe 

No. 

Plow i 
(1pm) 

I 

f 

length 

(m) 

J Segment 1 

I 

Segaent 2 

t 

1 « 

iDiameter ! 

! (mm) ; 

f t 

length 

(ml 

f 

iDiameter 
; (mm) 

f 

length 

(i^ 

1 

1200.00 

10000.0 

150 

1000.00 

- 

- 

2 

863.98 

1000.0 

150 

999 .03 

125 

0^97 

3 

199.59 

1000.0 

80 

877,46 

63 

122.54, 

4 

2,28 

1000.0 

20 

44.59 

15 

955.41 

5 

216.02 

1000.0 

80 

1000.00 

- 


6 

584.39 

1000.0 

125 

1000.00 


- 

7 

204.30 

1000.0 

80 

1000.00 

- 

- 

8 

1.87 

1000.0 

15 

1000.00 

- 

- 

9 

280.09 

1000.0 

100 

999.47 

80 

0.53 

10 

180.09 

1000.0 

80 

1000.00 

- 

- 

11 

11.15 

1000.0 

32 

1000.00 

- 

, 

12 

91 .24 

1000.0 

80 

. 14.01 

63 

985 t 99 

13 

8.76 

1000.0 

32 

1000.00 

- 

— 

14 

95.02 

1000.0 

63 

950.11 

50 

49«89 

15 

13.74 

1000.0 

32 

1000.00 

- 

— 

16 

100.00 

1000,0 

63 

1000.00 


— 

17 

13.74 

1000.0 

32 

1000.00 

- 

— 

18 

113.74 

1000.0 

63 

1000,00 


mm 

Optimum pumping 

head at 

node 1 = 

35 + 12, 

14 = 47 

.14 


Optimtam cost of the system = 3H11.34 Es/year 

Optimum cost of the initial flow patteami = 31852,17 Es/year 

No booster pump is required at pipe no. 6. 
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Table 23 

Optimum Design for Problem - 6 


p 

Pipe i 

No. ; 

1 

f 

f 

t 

Plow 

(1pm) 

i length 

1 (ml 

t 

I 

f 

s 

Segment 1 

; Segment 2 

I ■ 

J- 

IDiameter 

I* (mm) 1 

I 1 

length 

(m) 

""'"1 ' ' '1'"' ' ' ' ■■ 

Iliameter J length 
; (mm) ; (m) 

t f 

1 

1400.00 

1000.0 


150 

1000.00 


- 

2 

1084.71 

1000.0 


150 

1000.00 

- 


3 

199.18 

1000.0 


80 

684.31 

63 

351.69 

4 

2.03 

1000.0 . 


15 

1000.00 

- 

- 

5 

215.29 

1000.0 


80 

611.21 

63 

388.79 

6 

785.53 

1000.0 


125 

1000.00 

- 

- 

7 

312.06 

1000.0 


100 

1000.00 


- 

8 

1.58 

1000.0 


15 

1000.00 

- 

- 

9 

373.47 

1000.00 


100 

1000.00 

- 

- 

10 

273.47 

1000.0 


100 

977.19 

80 

22.81 

11 

5.05 

1000.0 


25 

9.69 

20 

990.31 

12 

178.52 

1000.0 


80 

1000.00 

- 

— 

13 

14.38 

1000.0 


32 

1000.00 

- 

— 

14 

208.59 

1000.0 


80 

1000.00 

— 

-- 

15 

12.89 

1000.0 


32 

1000.00 

- 

— 

16 

99.63 

1000.0 


63 

944.45 

50 

55.55 

17 

13.26 

1000.0 


32 

1000.00 

- 


18 

113.26 

1000.0 


63 

998.78 

50 

1,22 

19 

92.90 

1000.0 


80 

1000.00 

- 

— . 

20 

7.10 

1000.0 


32 

1000.00 

- 


21 

107.10 

1000.0 


80 

407.77 


592,23 


Optimum pumping head at node t = 35 * 24.40 = 59.40 Q 

Minimum cost of the system = 40174.68 Es/year 

Minimum cost for the initial flow pattern = 40813 *91 Rs/year 
Mo booster pump is required at pipe no. 6. 
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Table 24 

Optimtim Design for Problem - 7 


Pipe 

No. 


T 

t 

J length 

J (m) 

f 

t 

1 

1 

— r 

{ Segment 1 

I. 

# 

1 

Segment 2 

t 

, Diameter 

I (mm) 

I 

Length 

(ml 

f t 

Diameter 

I (mm) 1 

f ... . t 

Length 

(e1 

1 

1500.00 

1000.0 

150 

1000.00 


- 

- 

2 

1184.56 

1000.0 

150 

1000.00 


- 

- 

3 

308.49 

1000.0 

100 

499.32 


80 

500.68 

4 

2.80 

1000.0 

15 

1000.00 


- ^ , 


5 

215 .44 

1000.0 

80 

879.58 


63 

120.42 

6 

776.07 

1000.0 

125 

1000.00 


- 

- 

7 

300.06 

1000.0 

100 

962.91 


80 

37.09 

8 

2.01 

1000.0 

15 

1000.00 


- 

- 

9 

376.01 

1000.0 

100 

1000.00 


- 

- 

10 

276.01 

1000.0 

100 

1000.00 


- 

- 

11 

2.32 

1000.0 

20 

34.09 


15 

965.91 

12 

178.33 

1000.0 

80 

1000.00 


- 

■ ■- 

13 

15.34 

1000.0 

32 

1000.00 


- 

— 

14 

199.75 

1000.0 

80 

1000,00 


- 

- 

15 

11.32 

1000.0 

32 

837 .04 


25 

162.96 

16 

209.28 

1000.0 

80 

1000.00 


- 

- 

17 

12.64 

1000.0 

32 

1 000,00 


- 


18 

112.64 

1000.0 

63 

667.17 


50 

352.83 

19 

93.67 

1000.0 

80 

1000,00 


- 

. 

20 

6.32 

1000.0 

32 

1000.00 



■ 

21 

95.72 

1000.0 

63 

984 .96 


50 

15.0;34":;, 

22 

10.60 

1000,0 

32 

1000.00 


- ■' 


23 

110.60 

1000.0 

80 

23.21 


■63;,;' > 

976.79 


Optimum pumping head at node 1 = 35 * 26.89 = 61 .89 m 

Minimiam cost of the system = 43644.54 Ss/year 

Minimum cost for the initial flow distribution — 44399.74 Es /year 
No booster pump is required at pipe no. 6. 
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Table 25 

Optimtim Results for Problem - 8 


Pipe 

No. 

^ 

Plow ' 

(Ipm) ; 

f 

1 

len^h 

(m) 

— 1 f- 

; Segment 1 

f f 

Segment 2 

I 

IDiameter 

I (mm) 

1 

Len^h IDiameter J 

(m) ; (mm) I 

1 t 

Length 

Cm) 

1 

1600.00 

1000.0 

150 

1000.00 


' 

2 

1177.45 

1000.0 

150 

1000.00 

- 

- 

3 

307.09 

1000.0 

80 

1000.00 


w... 

4 

2.53 

1000.0 

15 

1000.00 

- 


5 

322.54 

1000.0 

100 

15.45 

80 

984.55 

6 

770.36 

1000.0 

125 

1000.00 

- 

- 

7 

292.13 

1000.0 

100 

556.24 

80 

443.76 

8 

9.43 

1000.0 

20 

947.68 

15 

52.32 

9 

378,23 

1000.0 . 

100 

1000.00 

- 

- 

10 

278.23 

1000,0 

100 

786.84 

80 

213.16 

11 

2,03 

1000.0 

15 

1000.00 


- 

12 

180.26 

1000.0 

80 

1000.00 

- 

- 

13 

1 .77 

1000.0 

20 

54.20 

15 

965 .80 

14 

199.53 

1000.0 

80 

1000.00 

- 

- , 

15 

11.50 

1000.0 

32 

1000.00 


— 

16 . 

200.19 

1000.0 

80 

1000.00 

- 

- 

17 

2.82 

1000.0 

20 

435.61 

15 

'564.38 

18 

220.01 ■ 

1000.0 

80 

1000.00 

- 

— 

19 

82.04 

1000.0 

63 

1000.00 



20 

17.96 

1000.0 

50 

1000.00 



21 

109.26 

1000.0 

80 

11 .61 

63 

^8.39 

22 

8.70 

1000.0 

32 

1000.00 

- 


23 

91.51 

1000.0 

63 

705.89 

50 

;;294.i;o;'; ■' 

24 

17.19 

1000.0 

40 

45.33 

32 

956':.6?' 

25 

117.19 

1000.0 

65 

1000.00 



Optimum pmping 

head at 

node 1 = 

35 + 32.51 

= 67 . 

51 m 


Minimum cost of the system = 46819.07 Rs/year 

Minimum cost for initial flow distribution = 47581.68 Rs/year 

No booster pump is required at pipe no. 9. 
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As mentioned earlier also, 15 mm dia pipe should 
not have been included in the list of available pipes but 
as it was included so the program in its attempt to minimize 
the objective function has provided 15 mm pipes for some 
minor pipe connections. 

The program gives the lengths of pipe segments of 
adjacent available diameters to be provided between any two 
nodes in odd figures. These should be rounded off to the 
nearest suitable values. Then a network solver should be 
run to check for the heads available at the nodes. This 
has not been done in the present case and the pipe lengths 
have been reported upto two decimal places (These have to 
be rounded off, if possible, to get the multiples of a 
preassigned value, say 1 meter). 

To check if the program worked well with respect 
to the pressure constraints and loop constraints, the head 
losses in pipes and pressures at nodes are computed for 
Problem 1 and are reported in Table 26. The pressure const- 
raints and loop constraints are found to be satisfied. 

The program ms executed for the 8 problems given 
in the preceding section. The execution time along with the 
size of problem is given in Table 27. The increase in 
Central Processing Unit (CPU) time with the increasing 
size of the problem is shown in Figure 10. 
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Table 26 

Head losses in Pipes and Available Heads at 
Nodes for Problem - 1 


T- 

1 

Pipe ; 

No., r 

t 

» 

Head loss (m) 

^ 

Segment 1 Segment 2 

f 

— 1 

i Node i 

*;no. ; 

f 1 

f f 

Head available 
(m) 

1 

16.417 

m 

1 

53.26 

2 

12.143 

- 

2 

36.844 

5 

5.426 

- 

3 

24.701 

4 

6.944 

- 

4 

26.219 

5 

5.886 

4.738 

5 

19.275 

6 

4.275 

- 

6 

18.953 

7 

3.953 

- 

7 

15.000 

8 

7.136 

0.131 




Table 27 

■Variation of Computer Time With Size of Problems Taken Up 


Problem 

t I J 

CPU time ko. of icFU time .'No. of 

No, of 

No. 

(seconds )4flow Jper Idecision 

{iteration {iteration {variables 

constraints 


{ {(seconds) ! 

, t f » 

f . f . t - . - - 



1 

3.08 

11 

0.28 

39 

12 

2 

19.31 

24 

0.804 

53 

17 

3 

9.71 

9 

1 .078 

62 

20 

4 

33.92 

16 

2.12 

76 

25 

5: ^ • 

54.23 

18 

3.013 

85 

28 

6 

68.13 

16 

4.258 

98 

32 

7 

121.54 

20 

6.077 

107 

35 

8 

81.59 

11 

7.414 

116 

38 
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6. comusiois AND SUGGESnOFS fOR PUTDRE WORK 
6 , 1 Conc lusions 

The program was found to work well for the 
problems framed in Section 4. The CPU time increases 
rapidly with the increasing size of the problem. The 
salient features of the program are as f ollowss 

(a) The program does not require a network solver. It 
incorporates the flow solution into the optimization 
procedure, v/ithout making any assumptions about the 
hydraulic solution of the network. 

(b) Operational decisions are included explicitly in 
the design process. 

(c) The method yields a design that is hydraulically 
feasible and is closer to being optimal than the 
one from which the search is started; this holds 
true even when the optimization procedure is 
terminated prematurely. 

(d) Algorithm can be used to satisfy exact flow require- 
ments along the pipes if required. 

Some weak points of the algorithm are; 

(a) Computer memory required is quite large as every 
pipe requires three decision variables in the IPG 
method . 

(b) Plows into and out of the reservoirs have to be fixed 
beforehand. This ensures proper operation of the 
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reservoir but does not include their capacity as a 
decision variable, 

^•2 Suggestions for further Work 

(1) The algorithm can be made more efficient by incorpo- 
rating some features as given belows 

a) The candidate list of diameters for any pipe is 
extended on the higher or lower side according 
to the requirement and the problem is reset and 
solved again. However if the number of diameter 
changes required is not significant compared to 
total number of decision variables, the post- 
optimality analysis can be used to get the 
solution using the solutions of the last problem. 

b) Similarly when flows in the pipes are changed, 
the post-optimality considerations can be 
employed to save computer time. 

c) A unidirectional search method can be suitably 
used to deteimine the step length while deciding 
about the flow change, 

(2) The program can be made to incorporate consideration 
of more than one demand pattern as proposed by 
Alperovits and Shamir^ by introducing 'artificial’ 
valves in each loop. 

(3) The program can be used for optimal expai^ion of an 
existing system or to decide about the operating 



policies of an existing system as outlined by 

0 

Alperovits and Shamir . 

There is no provision in the algorithm to reverse 
the flow directions in some pipes if required to 
minimize the cost. Hence some criterion should be 
developed to come out with a few suitable flow 
patterns and the one which gives the minimum cost 
should be taken as the optimum design. 

As computer memory required is quite large for large 
water distribution systems hence the application of 
sparce matrix techniques can save a considerable 
computer memory and time. 

More efficient' algorithms for network optimization 
may be tried (e.g. Out of Kilter Algorithm). 
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APPENDIX 1 

DERIVATION OP GRADIENT EXPRESSION 

Eor illustrative pixrpose the gradient expression 
given in general farm equation (19), is derived for the 
example problem shown in Pigure 11. For simplicity only 
the two loop paths are considered. It is assumed that 
additional., paths, which would be used to specify minimum 
node pressures, are not included, and therefore such paths 
are not considered in the derivation of the gradient expre- 
ssions. New notation is defined as follows s 
Hlj - =s head loss in the links between nodes i and j in the 

X J 

direction of the flow? 

= flow in the link between i and j? 
d X - change in flow in all the links in loop Aj 
d ^ = change in flow in all the links in loop 

y = head discontinuity at node 4 from loop A? 

S = head discontinuity at node 5 from loop B, 

The loop constraints from the linear program can 
be written ass 

Hl^^ — Hl^^ — Hl2'^ ~ ^ (50) 

- H1^5 * * Hlg.^ t = 0 (51) 

The dual variables for these constraints, and 
¥g, respectively, represent the change in cost, that would 
result if the head loss in the loops could change. If a 
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discontinuity in head loss is introduced , and can 
be used to determine whether the discontinuities in the 
ri^t hand sides of the constraints should be positive or 
negative to reduce the cost of the system. In the physical 
system, changes in the flows in the links are necessary to 
balance the head losses (i.e., to ’heal’ the discontinuity). 
Since flow and head losses are not independent, the neces- 
sary changes in flovns can be predicted. 

If the change in flow in each link is noted as 




- did 

(52) 

30-24 

= 

dfX 

(53) 



d^ - dfb 

(54 ) 

<3Q46 

:s 

d|3 

(55) 


= ■ 


(56) 



- di d- 

(57) 


= 

dp 

(58) 

The discontinuity expressions 

are written explicitly as 



HI 

24 = y 

(59) 

- HI 45 * * 

HI 

75 " ^ 

( 60 ) 

The total change in head loss 

around the loops due to 


t 

changes in flows should be 

- 

y and . 
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5 


3H1„3 

- 3Q23 '^'^24 


3045 


2H1.^ OHlf-ry ^Hl^c 

^<^46 ^ |^%7 ^ 


kP ^<^45 * ^0,, * 


Substituting for dQ. . from equations (52) to (58) 

^Hl,. 3 HI,, 

(dot - dft) - (- a=i) - (- dU) 

5*^45 “ ^%5 9^23 


3 ^ 2 ^ 
3%4 


dfl(. + dy = 0 


(63) 






= 0 


(64) 


All the information required to determine the change in 
cost with respect to a change in the flow in the paths is 
now available. 


^^j^ lGoSt 1 LVostJ + _£iiiO£Li SL2_ (6*51 

33 ^ Soi. 

3 (Cost) 3 (Cost) -2 y , 9 (Cost) 2^ ^ (66) 

" 3y ’ 9^ ' df> 

The terms 3(Cost)/9y and ^>(Cost)/o)$ are the 
dual variables for the loops, and the remaining terms can be 
found by picking out the coefficients of the differentials 
in equations (63) aii<3 (64). 
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ill . . 1!^ . '^“45 . . 9 
3^4 9^55 9^23 9^2>r 

ijL ^ ^ ^ ^^45 

Dp- a-^ 

ai) ^ 9 ^^45 

'3<A ~ 2%5 




(67) 

( 68 ) 

(69) 

(70) 


Substituting these partial derivatives and the dual 

variables into equations (65) and (66), and using equation 
3 HI. . 

((1 ) f or , we get 


(Cost) u /^45 ^35 + 5: 


^23 *^24 


Hl/ic: 


(71) 


2 (Cost 


HI.. HI,. HI.., HI. 


= VI 


• Q/ir B H/ic; 0/!A Ogy f-JY5 


(72) 


A generalization of the above expression leads to 
the gradient e quat ion given . previously as (19). 
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APPENDIX 2 
NOTATIONS 

Q ss Discharge from a pipe 

Cjj^ ss Hazen-Williains coefficient 

D = Diameter of pipe 

H ss Head loss between the ends oT a pipe 
L ss Length of pipe 

J = Hydraulic gradient 

a, b, c * Numerical coefficients 

i, j: = Nodes- 

m = Number of pipe diameter in the candidate 

list of diameters 


0 = Cost factor for the decision variable 

1 > Inflow 

X = Decision variable 




